Multiomics analysis reveals B. MO1 as a distinct Babesia species and provides insights into its evolution and virulence.

Babesiosis, caused by protozoan parasites of the genus Babesia , is an emerging tick-borne disease of significance for both human and animal health. Babesia parasites infect erythrocytes of vertebrate hosts where they develop and multiply rapidly to cause the pathological symptoms associated with the disease. The identification of various Babesia species underscores the ongoing risk of new zoonotic pathogens capable of infecting humans, a concern amplified by anthropogenic activities and environmental shifts impacting the distribution and transmission dynamics of parasites, their vectors, and reservoir hosts. One such species, Babesia MO1, previously implicated in severe cases of human babesiosis in the midwestern United States, was initially considered closely related to B. divergens , the predominant agent of human babesiosis in Europe. Yet, uncertainties persist regarding whether these pathogens represent distinct variants of the same species or are entirely separate species. We show that although both B. MO1 and B. divergens share similar genome sizes, comprising three nuclear chromosomes, one linear mitochondrial chromosome, and one circular apicoplast chromosome, major differences exist in terms of genomic sequence divergence, gene functions, transcription profiles, replication rates and susceptibility to antiparasitic drugs. Furthermore, both pathogens have evolved distinct classes of multigene families, crucial for their pathogenicity and adaptation to specific mammalian hosts. Leveraging genomic information for B. MO1, B. divergens , and other members of the Babesiidae family within Apicomplexa provides valuable insights into the evolution, diversity, and virulence of these parasites. This knowledge serves as a critical tool in preemptively addressing the emergence and rapid transmission of more virulent strains.

Rats exhibit age-related mosaic loss of chromosome Y.

Mosaic loss of the Y chromosome (LOY) is the most frequent chromosomal aberration in aging men and is strongly correlated with mortality and disease. To date, studies of LOY have only been performed in humans, and so it is unclear whether LOY is a natural consequence of our relatively long lifespan or due to exposure to human-specific external stressors. Here, we explored whether LOY could be detected in rats. We applied a locus-specific PCR and target sequencing approach that we used as a proxy to estimate LOY in 339 samples covering eleven tissues from young and old individuals. We detected LOY in four tissues of older rats. To confirm the results from the PCR screening, we re-sequenced 60 full genomes from old rats, which revealed that the Y chromosome is the sole chromosome with low copy numbers. Finally, our results suggest that LOY is associated with other structural aberrations on the Y chromosome and possibly linked to the mosaic loss of the X chromosome. This is the first report, to our knowledge, demonstrating that the patterns of LOY observed in aging men are also present in a rodent, and conclude that LOY may be a natural process in placental mammals.

High-throughput transcriptome sequencing and comparative analysis of Escherichia coli and Schizosaccharomyces pombe in respiratory and fermentative growth.

In spite of increased complexity in eukaryotes compared to prokaryotes, several basic metabolic and regulatory processes are conserved. Here we explored analogies in the eubacteria Escherichia coli and the unicellular fission yeast Schizosaccharomyces pombe transcriptomes under two carbon sources: 2% glucose; or a mix of 2% glycerol and 0.2% sodium acetate using the same growth media and growth phase. Overall, twelve RNA-seq libraries were constructed. A total of 593 and 860 genes were detected as differentially expressed for E. coli and S. pombe, respectively, with a log2 of the Fold Change ≥ 1 and False Discovery Rate ≤ 0.05. In aerobic glycolysis, most of the expressed genes were associated with cell proliferation in both organisms, including amino acid metabolism and glycolysis. In contrast in glycerol/acetate condition, genes related to flagellar assembly and membrane proteins were differentially expressed such as the general transcription factors fliA, flhD, flhC, and flagellum assembly genes were detected in E. coli, whereas in S. pombe genes for hexose transporters, integral membrane proteins, galactose metabolism, and ncRNAs related to cellular stress were overexpressed. In general, our study shows that a conserved "foraging behavior" response is observed in these eukaryotic and eubacterial organisms in gluconeogenic carbon sources.

Bacterial and Fungal Microbiome Profiling in Chilhuacle Negro Chili (Capsicum annuum L.) Associated With Fruit Rot Disease.

Chilhuacle negro chili (Capsicum annuum L.) is an ancient Mexican landrace that is deeply linked to the culinary heritage of the country. Because of the high profitability and uniqueness of this crop, the Universidad Autónoma del Estado de Morelos is exploring its production in controlled environments. In the crop cycles of 2018 to 2019, the production of chilhuacle negro plants was seriously affected by an unidentified pathogen causing fruit rot, which reduced its quality, yield, and market value. Therefore, the main objective of this work was to study and characterize the fruit microbiota, which could help reveal the causal agent of this disease. Using DNA metabarcoding coupled with Illumina and nanopore sequencing technologies, we collected and analyzed both healthy and infected chili fruit, along with greenhouse bioaerosols. We also explored the bacterial and fungal microbiota by using microbiological techniques to isolate some of the culturable bacterial and fungal species. Our results suggest that the seedborne fungus Alternaria alternata is activated during the maturation stage of chilhuacle negro fruit, triggering a microbiome imbalance, which may in turn enable the establishment of other opportunistic pathogenic fungi during fruit decay, such as Mucor sp. To our knowledge, this is the first study of the chilhuacle negro chili microbiome, which can shed some light on our understanding of one of the main diseases that affect this valuable crop.

Genomic Analysis of Early SARS-CoV-2 Variants Introduced in Mexico.

The coronavirus disease 2019 (COVID-19) pandemic has affected most countries in the world. Studying the evolution and transmission patterns in different countries is crucial to enabling implementation of effective strategies for disease control and prevention. In this work, we present the full genome sequence for 17 SARS-CoV-2 isolates corresponding to the earliest sampled cases in Mexico. Global and local phylogenomics, coupled with mutational analysis, consistently revealed that these viral sequences are distributed within 2 known lineages, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage A/G, containing mostly sequences from North America, and lineage B/S, containing mainly sequences from Europe. Based on the exposure history of the cases and on the phylogenomic analysis, we characterized 14 independent introduction events. Additionally, three cases with no travel history were identified. We found evidence that two of these cases represented local transmission cases occurring in Mexico during mid-March 2020, denoting the earliest events described for the country. Within this local transmission cluster, we also identified an H49Y amino acid change in the Spike protein. This mutation represents a homoplasy occurring independently through time and space and may function as a molecular marker to follow any further spread of these viral variants throughout the country. Our results provide a general picture of the SARS-CoV-2 variants introduced at the beginning of the outbreak in Mexico, setting the foundation for future surveillance efforts.IMPORTANCE Understanding the introduction, spread, and establishment of SARS-CoV-2 within distinct human populations as well as the evolution of the pandemics is crucial to implement effective control strategies. In this work, we report that the initial virus strains introduced in Mexico came from Europe and the United States and that the virus was circulating locally in the country as early as mid-March. We also found evidence for early local transmission of strains with a H49Y mutation in the Spike protein, which could be further used as a molecular marker to follow viral spread within the country and the region.

Comparative and functional genomics of the protozoan parasite Babesia divergens highlighting the invasion and egress processes.

Babesiosis is considered an emerging disease because its incidence has significantly increased in the last 30 years, providing evidence of the expanding range of this rare but potentially life-threatening zoonotic disease. Babesia divergens is a causative agent of babesiosis in humans and cattle in Europe. The recently sequenced genome of B. divergens revealed over 3,741 protein coding-genes and the 10.7-Mb high-quality draft become the first reference tool to study the genome structure of B. divergens. Now, by exploiting this sequence data and using new computational tools and assembly strategies, we have significantly improved the quality of the B. divergens genome. The new assembly shows better continuity and has a higher correspondence to B. bovis chromosomes. Moreover, we present a differential expression analysis using RNA sequencing of the two different stages of the asexual lifecycle of B. divergens: the free merozoite capable of invading erythrocytes and the intraerythrocytic parasite stage that remains within the erythrocyte until egress. Comparison of mRNA levels of both stages identified 1,441 differentially expressed genes. From these, around half were upregulated and the other half downregulated in the intraerythrocytic stage. Orthogonal validation by real-time quantitative reverse transcription PCR confirmed the differential expression. A moderately increased expression level of genes, putatively involved in the invasion and egress processes, were revealed in the intraerythrocytic stage compared with the free merozoite. On the basis of these results and in the absence of molecular models of invasion and egress for B. divergens, we have proposed the identified genes as putative molecular players in the invasion and egress processes. Our results contribute to an understanding of key parasitic strategies and pathogenesis and could be a valuable genomic resource to exploit for the design of diagnostic methods, drugs and vaccines to improve the control of babesiosis.

Evolutionary analysis of the Chikungunya virus epidemic in Mexico reveals intra-host mutational hotspots in the E1 protein.

The epidemic potential of Chikungunya virus (CHIKV) was recently made evident by its introduction and rapid expansion in the Caribbean and the Americas. We sought to gain a detailed understanding of the dynamics of the epidemic in Mexico, the country with the highest number of confirmed CHIKV cases in the Americas, and to characterise viral evolution at the population and intra-host level. Analysis of the spatiotemporal distribution of 2,739 diagnosed cases in Mexico from December 2014 to December 2015 showed a rapid nationwide expansion of the epidemic with focalisation in the South West of the country. We sequenced the envelope glycoprotein 1 gene (E1) from 25 patients using the Illumina MiSeq platform and report synonymous and non-synonymous consensus mutations. Bayesian phylogenetic analysis using 249 Asian lineage E1 sequences gave updated estimates of nucleotide substitution rates for E1 and time to most recent common ancestor of major lineages. The analysis indicates phylogenetically-related emergent Latin American clusters in South Western Mexico, Nicaragua and Honduras and transmission of American strains in the Pacific islands. Detailed analysis showed that intra-host changes in E1 mainly occurred in two variable regions (E1:189-220 and E1:349-358) in domains II and III, respectively, in residues involved in inter and intra-envelope spike interactions. At the population level, this study sheds light on the introduction and evolutionary dynamics of CHIKV in the Americas. At the intra-host level, this study identifies mutational hotspots of the E1 protein with implications for understanding the relationship between the CHIKV quasispecies, viral fitness and pathogenesis.

The C-terminal protoxin region of Bacillus thuringiensis Cry1Ab toxin has a functional role in binding to GPI-anchored receptors in the insect midgut.

Bacillus thuringiensis Cry toxins are used worldwide for controlling insects. Cry1Ab is produced as a 130-kDa protoxin that is activated by proteolytic removal of an inert 500 amino-acid-long C-terminal region, enabling the activated toxin to bind to insect midgut receptor proteins, leading to its membrane insertion and pore formation. It has been proposed that the C-terminal region is only involved in toxin crystallization, but its role in receptor binding is undefined. Here we show that the C-terminal region of Cry1Ab protoxin provides additional binding sites for alkaline phosphatase (ALP) and aminopeptidase N (APN) insect receptors. ELISA, ligand blot, surface plasmon resonance, and pulldown assays revealed that the Cry1Ab C-terminal region binds to both ALP and APN but not to cadherin. Thus, the C-terminal region provided a higher binding affinity of the protoxin to the gut membrane that correlated with higher toxicity of protoxin than activated toxin. Moreover, Cry1Ab domain II loop 2 or 3 mutations reduced binding of the protoxin to cadherin but not to ALP or APN, supporting the idea that protoxins have additional binding sites. These results imply that two different regions mediate the binding of Cry1Ab protoxin to membrane receptors, one located in domain II-III of the toxin and another in its C-terminal region, suggesting an active role of the C-terminal protoxin fragment in the mode of action of Cry toxins. These results suggest that future manipulations of the C-terminal protoxin region could alter the specificity and increase the toxicity of B. thuringiensis proteins.

Spodoptera frugiperda (J. E. Smith) Aminopeptidase N1 Is a Functional Receptor of the Bacillus thuringiensis Cry1Ca Toxin.

Bacillus thuringiensis Cry1Ca is toxic to different Spodoptera species. The aims of this work were to identify the Cry1Ca-binding proteins in S. frugiperda, to provide evidence on their participation in toxicity, and to identify the Cry1Ca amino acid residues involved in receptor binding. Pulldown assays using Spodoptera frugiperda brush border membrane vesicles (BBMV) identified aminopeptidase N (APN), APN1, and APN2 isoforms as Cry1Ca-binding proteins. Cry1Ca alanine substitutions in all residues of domain III ß16 were characterized. Two ß16 nontoxic mutants (V505A and S506A) showed a correlative defect on binding to the recombinant S. frugiperda APN1 (SfAPN1). Finally, silencing the expression of APN1 transcript, by double-stranded RNA (dsRNA) feeding, showed that silenced larvae are more tolerant of the Cry1Ca toxin, which induced less than 40% mortality in silenced larvae whereas nonsilenced larvae had 100% mortality. Overall, our results show that Cry1Ca relies on APN1 binding through domain III ß16 to impart toxicity to S. frugiperdaIMPORTANCEBacillus thuringiensis Cry toxins rely on receptor binding to exert toxicity. Cry1Ca is toxic to different populations of S. frugiperda, a major corn pest in America. Nevertheless, the S. frugiperda midgut proteins that are involved in Cry1Ca toxicity have not been identified. Here we identified aminopeptidase N1 (APN1) as a functional receptor of Cry1Ca. Moreover, we showed that Cry1Ca domain III ß16 is involved in APN1 binding. These results give insights on potential target sites for improving Cry1Ca toxicity to S. frugiperda.

Genome Sequence of Torulaspora delbrueckii NRRL Y-50541, Isolated from Mezcal Fermentation.

Torulaspora delbrueckii presents metabolic features interesting for biotechnological applications (in the dairy and wine industries). Recently, the T. delbrueckii CBS 1146 genome, which has been maintained under laboratory conditions since 1970, was published. Thus, a genome of a new mezcal yeast was sequenced and characterized and showed genetic differences and a higher genome assembly quality, offering a better reference genome.

High-Quality Draft Genome Sequence of Candida apicola NRRL Y-50540.

Candida apicola, a highly osmotolerant ascomycetes yeast, produces sophorolipids (biosurfactants), membrane fatty acids, and enzymes of biotechnological interest. The genome obtained has a high-quality draft for this species and can be used as a reference to perform further analyses, such as differential gene expression in yeast from Candida genera.

High-Quality Draft Genome Sequence of Babesia divergens, the Etiological Agent of Cattle and Human Babesiosis.

Babesia divergens causes significant morbidity and mortality in cattle and splenectomized or immunocompromised individuals. Here, we present a 10.7-Mb high-quality draft genome of this parasite close to chromosome resolution that will enable comparative genome analyses and synteny studies among related parasites.

RegulonDB version 7.0: transcriptional regulation of Escherichia coli K-12 integrated within genetic sensory response units (Gensor Units).

RegulonDB (http://regulondb.ccg.unam.mx/) is the primary reference database of the best-known regulatory network of any free-living organism, that of Escherichia coli K-12. The major conceptual change since 3 years ago is an expanded biological context so that transcriptional regulation is now part of a unit that initiates with the signal and continues with the signal transduction to the core of regulation, modifying expression of the affected target genes responsible for the response. We call these genetic sensory response units, or Gensor Units. We have initiated their high-level curation, with graphic maps and superreactions with links to other databases. Additional connectivity uses expandable submaps. RegulonDB has summaries for every transcription factor (TF) and TF-binding sites with internal symmetry. Several DNA-binding motifs and their sizes have been redefined and relocated. In addition to data from the literature, we have incorporated our own information on transcription start sites (TSSs) and transcriptional units (TUs), obtained by using high-throughput whole-genome sequencing technologies. A new portable drawing tool for genomic features is also now available, as well as new ways to download the data, including web services, files for several relational database manager systems and text files including BioPAX format.

Genome-wide identification of transcription start sites, promoters and transcription factor binding sites in E. coli.

Despite almost 40 years of molecular genetics research in Escherichia coli a major fraction of its Transcription Start Sites (TSSs) are still unknown, limiting therefore our understanding of the regulatory circuits that control gene expression in this model organism. RegulonDB (http://regulondb.ccg.unam.mx/) is aimed at integrating the genetic regulatory network of E. coli K12 as an entirely bioinformatic project up till now. In this work, we extended its aims by generating experimental data at a genome scale on TSSs, promoters and regulatory regions. We implemented a modified 5' RACE protocol and an unbiased High Throughput Pyrosequencing Strategy (HTPS) that allowed us to map more than 1700 TSSs with high precision. From this collection, about 230 corresponded to previously reported TSSs, which helped us to benchmark both our methodologies and the accuracy of the previous mapping experiments. The other ca 1500 TSSs mapped belong to about 1000 different genes, many of them with no assigned function. We identified promoter sequences and type of sigma factors that control the expression of about 80% of these genes. As expected, the housekeeping sigma(70) was the most common type of promoter, followed by sigma(38). The majority of the putative TSSs were located between 20 to 40 nucleotides from the translational start site. Putative regulatory binding sites for transcription factors were detected upstream of many TSSs. For a few transcripts, riboswitches and small RNAs were found. Several genes also had additional TSSs within the coding region. Unexpectedly, the HTPS experiments revealed extensive antisense transcription, probably for regulatory functions. The new information in RegulonDB, now with more than 2400 experimentally determined TSSs, strengthens the accuracy of promoter prediction, operon structure, and regulatory networks and provides valuable new information that will facilitate the understanding from a global perspective the complex and intricate regulatory network that operates in E. coli.

Sensitive genome-wide screen for low secondary enzymatic activities: the YjbQ family shows thiamin phosphate synthase activity.

Contemporary enzymes are highly efficient and selective catalysts. However, due to the intrinsically very reactive nature of active sites, gratuitous secondary reactions are practically unavoidable. Consequently, even the smallest cell, with its limited enzymatic repertoire, has the potential to carry out numerous additional, very likely inefficient, secondary reactions. If selectively advantageous, secondary reactions could be the basis for the evolution of new fully functional enzymes. Here, we investigated if Escherichia coli has cryptic enzymatic activities related to thiamin biosynthesis. We selected this pathway because this vitamin is essential, but the cell's requirements are very small. Therefore, enzymes with very low activity could complement the auxotrophy of strains deleted of some bona fide thiamin biosynthetic genes. By overexpressing the E. coli's protein repertoire, we selected yjbQ, a gene that complemented a strain deleted of the thiamin phosphate synthase (TPS)-coding gene thiE. In vitro studies confirmed TPS activity, and by directed evolution experiments, this activity was enhanced. Structurally oriented mutagenesis allowed us to identify the putative active site. Remote orthologs of YjbQ from Thermotoga, Sulfolobus, and Pyrococcus were cloned and also showed thiamin auxotrophy complementation, indicating that the cryptic TPS activity is a property of this protein family. Interestingly, the thiE- and yjbQ-coded TPSs are analog enzymes with no structural similarity, reflecting distinct evolutionary origin. These results support the hypothesis that the enzymatic repertoire of a cell such as E. coli has the potential to perform vast amounts of alternative reactions, which could be exploited to evolve novel or more efficient catalysts.

Specific epitopes of domains II and III of Bacillus thuringiensis Cry1Ab toxin involved in the sequential interaction with cadherin and aminopeptidase-N receptors in Manduca sexta.

The Bacillus thuringiensis Cry toxins are specific to different insects. In Manduca sexta cadherin (Bt-R1) and aminopeptidase-N (APN) proteins are recognized as Cry1A receptors. Previous work showed that Cry1Ab binds to Bt-R1 promoting the formation of a pre-pore oligomer that binds to APN leading to membrane insertion. In this work we characterized the binding epitopes involved in the sequential interaction of Cry1Ab with Bt-R1 and APN. A Cry1Ab immune M13 phage repertoire was constructed using antibody gene transcripts of bone marrow or spleen from a rabbit immunized with Cry1Ab. We identified antibodies that recognize domain II loop 3 (scFvL3-3) or beta16-beta22 (scFvM22) in domain III. Enzyme-linked immunosorbent assay and toxin overlay binding competition assays in the presence of scFvL3-3, scFvM22, or synthetic peptides showed that domain II loop 3 is an important epitope for interaction with Bt-R1 receptor, whereas domain III beta16 is involved in the interaction with APN. Both scFvL3-3 and scFvM22 lowered the toxicity of Cry1Ab to M. sexta larvae indicating that interaction with both receptors is important for in vivo toxicity. scFvL3-3 and anti-loop2 scFv (scFv73) promoted the formation of the pre-pore oligomer in contrast to scFvM22. In addition, scFvL3-3 and scFv73 preferentially recognized the monomeric toxin rather than the pre-pore suggesting a conformational change in domain II loops upon oligomerization. These results indicate for the first time that both receptor molecules participate in Cry1Ab toxin action in vivo: first the monomeric toxin binds to Bt-R1 through loops 2 and 3 of domain II promoting the formation of the pre-pore inducing some structural changes, then the pre-pore interacts with APN through beta-16 of domain III promoting membrane insertion and cell death.

Nitric oxide in B6 mouse and nitric oxide-sensitive soluble guanylate cyclase in cat modulate acetylcholine release in pontine reticular formation.

ACh regulates arousal, and the present study was designed to provide insight into the neurochemical mechanisms modulating ACh release in the pontine reticular formation. Nitric oxide (NO)-releasing beads microinjected into the pontine reticular formation of C57BL/6J (B6) mice significantly (P < 0.0001) increased ACh release. Microdialysis delivery of the NO donor N-ethyl-2-(1-ethyl-2-hydroxy-2-nitrosohydrazino)-ethanamine (NOC-12) to the mouse pontine reticular formation also caused a concentration-dependent increase in ACh release (P < 0.001). These are the first neurochemical data showing that ACh release in the pontine reticular formation of the B6 mouse is modulated by NO. The signal transduction cascade through which NO modulates ACh release in the pontine reticular formation has not previously been characterized. Therefore, an additional series of studies quantified the effects of a soluble guanylate cyclase (sGC) inhibitor, 1H-[1,2,4]oxadiazolo-[4,3-a]quinoxalin-1-one (ODQ), on ACh release in the cat medial pontine reticular formation. During naturally occurring states of sleep and wakefulness, but not anesthesia, ODQ caused a significant (P < 0.001) decrease in ACh release. These results show for the first time that NO modulates ACh in the medial pontine reticular formation of the cat via an NO-sensitive sGC signal transduction cascade. Isoflurane and halothane anesthesia have been shown to decrease ACh release in the medial pontine reticular formation. The finding that ODQ did not alter ACh release during isoflurane or halothane anesthesia demonstrates that these anesthetics disrupt the NO-sensitive sGC-cGMP pathway. Considered together, results from the mouse and cat indicate that NO modulates ACh release in arousal-promoting regions of the pontine reticular formation via an NO-sensitive sGC-cGMP pathway.